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Abstract 

We discuss various physics aspects of neutrino oscillation with non-standard interactions (NSI) . 
We formulate a perturbative framework by taking Am2i/Am3]^, S13, and the NSI elements 8^0 
(a, (5 = e, /i, r) as small expansion parameters of the same order e. Within the e perturbation 
theory we obtain the S matrix elements and the neutrino oscillation probability formula to second 
order (third order in related channels) in e. The formula allows us to estimate size of the 
contribution of any particular NSI element e^^ to the oscillation probability in arbitrary channels, 
and gives a global bird-eye view of the neutrino oscillation phenomena with NSI. Based on the 
second-order formula we discuss how all the conventional lepton mixing as well as NSI parameters 
can be determined. Our results shows that while ^13, (5, and the NSI elements in v^, sector can in 
principle be determined, complete measurement of the NSI parameters in the — Uj- sector is not 
possible by the rate only analysis. The discussion for parameter determination and the analysis 
based on the matter perturbation theory indicate that the parameter degeneracy prevails with 
the NSI parameters. In addition, a new solar-atmospheric variable exchange degeneracy is found. 
Some general properties of neutrino oscillation with and without NSI are also illuminated. 
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I. INTRODUCTION 



Neutrino masses and lepton flavor mixing [1] discovered by the atmospheric [2], the 
solar and the reactor neutrino j3] experiments constitute still the uniques evidence for 
physics beyond the Standard Model. A possible next step would be a discovery of neutrino 
interactions outside the standard electroweak theory. Based on expectation of new physics 
at TeV scale such non-standard interactions (NSI) with matter possessed by neutrinos are 
proposed and extensively discussed [l, 0, 0, M, M, ID]. The experimental constraints on NSI 
are summarized in See also [l^ . 

Recognition of structure of neutrino masses and lepton flavor mixing, at least up to 
now, relies on neutrino flavor transformation [H, H, 13, 14], which we generically refer as 



neutrino oscillation in this paper. Quite naturally, there have been numerous theoretical 
analyses to understand the structure of the phenomena. In the context of long-baseline 
neutrino experiments, an exact expression of the oscillation probability is derived under the 



constant matter density approximation [15j. To understand physics of neutrino oscillation, 
however, it is often more illuminating to have suitable approximation schemes. In the 
latter category, various perturbative formulations of three-flavor neutrino oscillation have 
been developed and proven to be quite useful in particular in the context of long-basehne 
accelerator and reactor experiments. They include one-mass scale dominance approximation 
in vacuum [l6|, short- distance expansion in matter [13], matter perturbation theory [isl. [l9| . 
and perturbation theory with the small expansion parameters Am^i/Amli and ^13 [20(] 
(that are taken as of order e) which we call the e perturbation theory in this paper. See, 
for example, [2l|, 22, [2^, 24] for subsequent development of perturbation theory of neutrino 
oscillation. 

When the effects of NSI are included, however, theoretical analysis of the system of neu- 
trino flavor transformation does not appears to achieve the same level of completeness as 
that only with standard interactions (SI). Perturbative formulas of the oscillation proba- 



bilities with NSI have been derived under various assumptions [25|, l26|, |27|, l28|, l29|, l30|, l31 

Even some exact formulas are known jsij. However, one cannot answer the questions such 
as: How large is the effects of e^r in the oscillation probability P(z/e — f/^)? (See below for 
definition of NSI elements eap-) How large is the effects of in the oscillation probability 
P(z/e z/g)? Which set of measurement is sufficient to determine all the NSI elements? 
References of neutrino oscillation and the sensitivity analyses with NSI are too numerous 
to quote here and may be found in bibliographies in the existing literatures, for example, in 



It is the purpose of this paper to fill the gap between understanding of neutrino oscillations 
with and without NSI. We try to do it by formulating the similar e perturbation theory as 
in 20] but with including effects of NSI by assuming that NSI elements are of order ~ e. 
We derive the perturbative formula of the oscillation probability similar to the one in (20| . 
which we call, respectively, the NSI and the SI second-order formulas in this paper. The 
approximate formula will allow us to have a bird-eye view of the neutrino oscillations with 
NSI, and will enable us to answer the above questions. 

The other limitation that is present in some foregoing analyses, which we want to over- 
come, is the assumption of single (or, a few) eais dominance. Upon identification or getting 
hint for possible NSI interactions it will become possible to express Ea/B in propagation as 
a function of couplings involved in the higher dimensional operators. When this situation 
comes it is likely that all (or at least most of) the NSI elements exist in the Hamilto- 
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nian with comparable magnitudes. Therefore, the theoretical machinery we prepare for the 
analysis must include all the NSI elements at the same time.^ 

More about necessity and usefulness of the S matrix and the NSI second-order formula 
of all oscillation channels and with all NSI elements included; If we are to include the effects 
of NSI in production and detection processes it is necessary to sum up all the oscillation 
channels that can contribute. Hence, the formulas of all channels are necessary. In a previous 



paper it was uncovered that the so called ^13— NSI confusion [26|, |27j can be resolved by a 
two-detector setting in neutrino factory experiments lH. Keeping the terms with the solar 
Am^i is shown to be crucial for resolving the confusion, and hence a full second-order formula 
is useful. In fact, the NSI second-order formula is surprisingly simple in its form, keeping 
the form of the original SI one with generalized variables, and the structure is even more 
transparent than those with first-order approximation of NSI. 

With the NSI second order formula, we are able to discuss, for the first time, a strategy 
for simultaneous complete determination of the SI and NSI parameters. Through the course 
of discussions we indicate that, as in the system without NSI, the parameter degeneracy 



33l . |3J, |35| exist in systems with NSI, but in a new form which involve both the SI and the 
NSI parameters. See Sees. IVIII and lVIlH Moreover, we will uncover a new type of degeneracy, 
the one exchanging the generalized solar and atmospheric variables in Sec. IVII G[ 

Finally, we should mention about what will not be achieved in this paper even within 
the context of theoretical analysis. First of all, our perturbative formulation relies on the 
particular assumption on relative magnitudes of SI and NSI parameters, and we cannot say 
many for cases in which our assumptions are not valid. We discuss the effects of NSI while 
neutrinos propagate in matter, and its effects in production and detection of neutrinos are 
ignored. Therefore, this paper must be regarded as merely the first step toward complete 
treatment of neutrino oscillation with NSI. 



II. PHYSICS SUMMARY 



Because this paper has been developed into a long one, unfortunately, we think it con- 
venient for readers, in particular experimentalists, to summarize the physics outputs of the 
perturbative treatment of neutrino oscillation with NSI. We highly recommend the readers 
to read this section first. 



A. New result in the standard three flavor mixing 

Though this paper aims at uncovering structure of neutrino oscillation with NSI, we 
have observed a new features of standard neutrino oscillation without NSI in Sec. IV C\ the 
property we call the "matter hesitation" . It states that in our perturbative framework the 
matter effect comes in into the oscillation probability only at the second order in the small 
expansion parameter e in all the channels of neutrino oscillation.^ It is a highly nontrivial 



^ The similar comments also apply to the procedure by which the current constraints on NSI is derived 
(for example in 11]) where the constraints are derived under the assumption of presence of a particular 
NSI element in each time, the point carefully mentioned by the authors themselves. 

^ Though this property must be known in the community as the results of perturbative calculation, it 
appears to us that it did not receive enough attention so far. 
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feature because we treat the matter effect as of order unity. Tlie "matter liesitation" explains 
why it is so difficult to have a sufficiently large matter effect, e.g., to resolve the mass 
hierarchy, in many long-baseline neutrino oscillation experiments. It also has implications 
to neutrino oscillation with NSI as will be discussed in Sec. IV D[ 

B. Guide for experimentalists; importance of various NSI elements in each channel 

Experimentalists who want to hunt NSI in neutrino propagation may ask the following 
questions: 

• We want to uncover the effect of e^t (or ^e/x) • What is the neutrino oscillation channel 
do you recommend to use for this purpose? 

• We plan to detect the effect of e^r- Which set of measurement do we need to prepare? 

• We seek a complete determination of all the SI and the NSI parameters. What would 
be the global strategy to adopt? 

With the oscillation probability formulas given in Sec. IVII we will try to answer these ques- 
tions. Though we can offer only a partial answer to the last question above we can certainly 
give the answer to the first two questions within the framework of perturbation theory we 
use. In Table [T] the relative importance of the effects of each element Sap of NSI are tabu- 
lated as order of a small parameter e that they first appear in each oscillation channel. We 
presume e ~ 10~^. Thus, our answer to the above questions based on the assumption that 
only the terms up to second order in e are relevant would be (in order): 

• The neutrino oscillation channels in which only e^t and ^e/x come in and the other 
elements do not are the Ve related ones, z/g ^ z/g, t'e — ^^i-, and z/g — > Vr- Obviously, 
the latter two appearance channels would be more interesting experimentally. One 
can in principle determine them simultaneously with 6'i3 and 5 by rate only analysis. 

• Do measurement at the z/^ — z/^ channel to determine e^^. Adding z/^ — >■ channel 
does not help. Effects of the NSI element are relatively large because they are first 
order in e, but the spectrum information is crucial to utilize this feature and to separate 
its effects from those of — e-rr- If an extreme precision is required you might want 
to supplement the measurement by the and Vr appearance measurement above. 

• We will show that, in fact, there is a difficulty in complete determination of all the 
NSI and SI parameters by the rate only analysis. The trouble occurs in the — Vr 
sector. Even though we are allowed to assume perfect measurement of all the channels 
including the one with Vr beam (which, of course, would not be practical), one of the 
three unknowns, e^^ — Srr and e^r including its phase, cannot be determined if we 
rely on the rate only analysis. See Sec. IVIII for more details. Clearly, the spectrum 
information is the key to the potential of being able to determine all the SI and the 
NSI parameters, which should be taken into account in considering future facilities 
which search for NSI. 
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TABLE I: Presented are the order in e (~ 10"^) at which each type of = e,fi,T) and a 

dependence (a is Wolfenstein's matter effect coefficient 0]) starts to come in into the expression 
of the oscillation probability in e perturbation theory!^ The last column is for the a dependence 
in the standard oscillation without NSI. The order of e indicated in parentheses implies the one 
for the maximal O23 in which cancellation takes place in the leading order. See the text for the 
definition of e perturbation theory and for more details. 
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To second order in e the sensitivity to e^^ and Stt is through the form e^^ — Ett, and hence no sensitivity 
to the individual e's. Generally, the diagonal e's appear in a form of difference in the oscillation probabilities 
because an over-all phase is an unobservable. 



With regard to the second point above, some remarks are in order; Usually, disappearance 
channels are disadvantageous in looking for a small effect such as 6'i3, because one has to 
make the statistical error smaller than the effect one wants to detect. In this respect, the 
NSI search in the z/^ — z/^ sector is promising because it is the first order effect in e. In fact, 
rather high sensitivities for determining and — Err observed in atmospheric [36] and 
future accelerator |[37|] neutrino analyses are benefited by this feature. 

We must warn the readers that experimental observable will be affected by NSI effects in 
production and detection of neutrinos. Therefore, our comments in this subsection assumes 
that they are well under control and shown to be smaller than the NSI effects in propagation 
by near detector measurement with an extreme precision. It should be also emphasized 
that some of our comments rely on the second-order perturbative formula of the oscillation 
probability. 



C. Some interesting or peculiar features of neutrino oscillation with NSI 

We list here some interesting features of neutrino oscillation with NSI which will be fully 
discussed in the following sections in this paper. Some of them are either unexpected, or 
might be showed up in previous analyses but without particular attention. A few points in 
them requires further investigation for full understanding. 

• One of the most significant feature in Table. 1 is that See appears only at third order 
in e in all oscillation channels. It will be shown in Sec. IV Dl that this feature can be 
explained as a consequence of the matter hesitation mentioned earlier. 

• It is interesting to observe from Table [T] that Wolfenstein's matter effect coefficient 
a in the oscillation probability, shows up in first (second) order in e in system with 
(without) NSI, which makes effects of matter density uncertainty larger in system with 
NSI. It occurs in the z/^ — I'r sector, and is easily understood as a consequence of "tree 
level" transition by the NSI element. 

• The results in the last column in Table [T] indicates that sensitivity to — e^r will 
depend upon if 6*23 is maximal or not. This feature is clearly seen, e.g., in 33]. Analysis 
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to resolve the ^23 octant degeneracy similar to the one proposed for cases without NSI 
would be required for correct estimation of the sensitivity to NSI. 



III. INTRODUCING THE EFFECTS OF NSI IN NEUTRINO PRODUCTION, 
PROPAGATION AND DETECTION PROCESSES 

We consider NSI involving neutrinos of the type 

= -2v^5ijG^(z7,7M^L^/3) iJYPf), (1) 

where Gp is the Fermi constant, and / stands for the index running over fermion species in 
the earth, f = e,u, d, where P stands for a projection operator and is either = ^{1 — 75) 
or Pr = |(1 + 75). The current constraints on e parameters are summarized in [11']. 

Upon introduction of the NSI as in ([1]) it affects neutrino production, detection as well 
as propagation in matter [1, H^, 27, 2^. Therefore, we have to analyze the following "grand 



transition amplitude" from a parent 11 particles (which needs not to be pions) to the par- 
ticular detection particle N (which needs not be nucleons):^ 

T{Eu,En) = ^ / dE^aD{En, Eua)S{i'a — ^ ^13', Ei,a)R{E^i3, En) (2) 

where the sum over a and (3 must be taken only if they are amenable to be produced by the 
decay, or to undergo the reaction. Here, we have assumed the particular decay process to 
produce neutrinos as 11 ^ z/q, + Xa with decay amplitude D{Eyi, Eya) with the energies En 
and Eua of parent and daughter particles, and the particular reaction + Ptg — *■ Np + Yp 
with reaction amplitude R{E^p, En) which produce Np particle with energy E^- Here, Xa 
and Yp are meant to be some inclusive collections of particles and Ptg denotes the target 
particle. S'(z/q — > up) denotes the neutrino oscillation amplitude of the channel z/^ — > up. 
The observable quantity is of course \T{Eyi, En)\'^- 

We assume that the coupling constant Sap possessed by NSI is small, ~ (^'M^^ where 
Mjyp is a new physics scale, so that we can organize perturbative treatment of the effects of 
NSI. Eap can be as small as 10~^ (10~^) for M^p = 1(10) TeV, and is even smaller if higher 
dimension operators (higher than six) are required. We assume that all the e^p have similar 
order of magnitudes and denote the small number collectively as e. Under these assumptions 
we expect that the decay and the detection functions, and the oscillation probabilities can 
be expanded as 

D{Eu,E,^) = D{o) + DWe + D(V + ... 
S{u^ ^ up; E,^) = 5(°)(z/« ^ vp) + 5(i)(^a ^ + S^^\v^ ^ vp)e^ + ... 

R{E,^, En) = R^'^ + R^'h + R^^^e^ + ... (3) 

where we have suppressed the kinematical dependences in quantities in the right-hand-side 
of (13]). The first terms in ([3]) are the one without NSI. Now, because of the smallness of 



^ One can talk about momentum reconstructed detected neutrinos instead of detected positrons, for ex- 
ample, but the reconstruction process must involve the effects of NSI. The expression in ([2]) is just to 
symbolically indicate this point. 
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e ~ 10 ^ (or smaller) we take the attitude that keeping terms up to second order in e must 
be good enough to discuss the effects of NSI and eventually to estimate the sensitivity to 

Unfortunately, even with the perturbative treatment this is a highly complicated system 
to analyze its full structure. It is possible that the types of NSI that contribute to production 
and detection processes are more numerous than the ones in the propagation process [so]. If 
this occurs the effects of NSI into production and detection processes could be qualitatively 
different from those in propagation. Therefore, the effects of NSI come into the decay and 
the reaction amplitudes generally in a model-dependent fashion, so that the flavor {a, (3) 
dependence of NSI effects are also model-dependent. Also they do so in an energy dependent 
way so that integration over neutrino energy in ([2]) is required for the full analysis. For an 
explicit example of how NSI enter into the decay and the reaction amplitudes as well as 
to the neutrino propagation in matter in concrete models, see for example the "unitarity 



violation" approach developed in [43 



In this paper, therefore, we confine ourselves to analysis of the structure of neutrino 
propagation with NSI, namely the terms with no effects of NSI in the decay and the reaction 
amplitudes in ([2]). This is a particularly simple system (relatively speaking with the full one) 
in the sense that no unitarity violation comes in because it deals with propagation of three 
light neutrinos. Furthermore, it has no explicit model dependence once the effects of NSI is 
parametrized in the familiar way. See the Hamiltonian in ([5]). We should emphasize that 
limitation of our scope to the problem of neutrino propagation, in fact, allows us to dig out 
structure of neutrino oscillation with NSI in a transparent manner. Therefore, we think it 
a meaningful first step. 

Our analysis can become the whole story provided that extremely stringent bounds on 
NSI effects in decay as well as detection reactions are placed by front detector measurement 
in future experiments. Otherwise, it covers only a leading (zeroth) order terms in NSI effect 
in decay and detection. When the first order corrections to them are taken into account 
what is needed is to compute the oscillation amplitude up to first order in e to obtain the 
observable to order e^. Hence, we present the results of S matrix elements in Appendix \N\ 
not only the expression of the oscillation probabilities, for future use. 



IV. GENERAL PROPERTIES OF NEUTRINO OSCILLATION WITH AND 
WITHOUT NSI 

Now, we analyze the structure and the properties of neutrino propagation in matter with 
NSI. We, however, sometimes go back to the system without NSI whenever it is illuminating. 
The results obtained in this section are exact, that is, they are valid without recourse to 
perturbation theory we will formulate in the next section. To discuss effects of NSI on 
neutrino propagation it is customary to introduce the e parameters, which are defined as 
eai3 = S/p^^a^) where Uf {rie) denotes the /-type fermion (electron) number density 
along the neutrino trajectory in the earth. Then, the neutrino evolution equation can be 



^ As far as the appearance channels Ve — > i^fj. and i^e — > i^t are concerned the oscillation amplitudes start 
from first order in e, as we will see below. Therefore, only the first order corrections to D and R are 
relevant for the observable to order e^. 
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written in flavor basis as 



(4) 



In the standard three-flavor neutrino scheme, Hamiltonian including NSI is given by 



H = — (U 
2E 














' 1 
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+ a{x) 
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+ a(x) 



-i<f>e 



\^ flT I ^ 



(5) 



where Am^j = — mf, and a{x) = 2\/2GpNf,{x)E is the coeflicient which is related to the 
index of refraction of neutrinos in medium of electron number density Ng{x) [5|, where Gp is 
the Fermi constant and E is the neutrino energy. The flrst two terms in ([5]) are the Standard 
Model interactions, whereas the last term denotes the non-standard neutrino interactions 
with matter. U denotes the flavor mixing matrix, the Maki-Nakagawa-Sakata (MNS) matrix 
[l|, in the lepton sector. In its standard form [ij] it reads 



U = U23UnUi2 





—S23 C23 



ci3 sise 
1 

JS Q 



-iS 



C12 S12 
-S12 C12 
1 



(6) 



-sise- u Ci3 

where S stands for the leptonic Kobayashi-Maskawa (KM) phase j45^. 

Most of the formulas in this and the next sections (Sees. IIVI and |V]) can be written in 
forms valid for arbitrary matter density proflles if the adiabatic approximation holds. We, 
however, present the ones derived under the constant matter density approximation because 
it makes the equations simpler, in particular, the perturbative formulas for the oscillation 
probabilities in Sec. IVII Unlike the case of the MSW solar neutrino solutions in which 
the matter density variation is the key to the problem, the constant density approximation 
to Ne{x) in long-baseline experiments should serve as a reasonable flrst approximation. 

The 5* matrix describes possible flavor changes after traversing a distance L, 

Va{L) = Saf3l^f3{0), (7) 

and the oscillation probability is given by 

P{iyf,^u^;L) = \Saf,\^. (8) 
If the neutrino evolution is governed by the Schrodinger equation (jlj), S matrix is given as 



S = Texp 



—t I dxH{x) 





(9) 



where T symbol indicates the "time ordering" (in fact "space ordering" here). The right- 
hand-side of dH]) may be written as e~'^^^ for the case of constant matter density. For 
notational convenience, we denote the S matrix elements as 



S 



acq 

^ee *-^e/x ^er 
C q C 

Sre ^T/i S-j-f 



(10) 
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The primary purpose of this paper is to discuss the properties of neutrino oscillation 
in the standard three flavor system with NSI. But, we recollect the properties of neutrino 
oscillation without NSI whenever necessary, and treat both systems simultaneously or go 
back and forth between them to make our discussion transparent. By this way the properties 
of the neutrino oscillations can be better illuminated. 



A. Relations between neutrino oscillation amplitudes without NSI 



If NSI, the third term in (Q, is absent the matter term (the second term in ([5])) has a 
symmetry; It is invariant under U23 rotation which act on z/^ — subspace. Due to this 
symmetry the Hamiltonian can be conveniently written in the form 

H = U23HUI, (11) 

and hence the S matrix can be written as 

S{L) = U23SiL)Ul, 

as noticed in 2^ where S{L) = Texp —i dxH{x) . 
do not contain 923. 

If we denote S{L) matrix elements in a form analogous to the one in (fTOj) S matrix can 
be written as 

See 

C23Sij.e + S23ST 



(12) 

The point here is that H and S{L) 



C23'S'e^ + S23S er 
S2^SrT + C23S23{SfiT + Srfi) 



-S23Se^l + C23S er 

■523'S'r/x + C23S23{Stt — Sf^f^) 



. ~S23Sfie + C23'S'-re ^238-^^ — ■523'S'^t + C23'523('S'tt ~ 'S'^m) -^23^1^1^ + ^235*,-^ — C23S23(>S'^r + St^^ 



It should be noticed that See = See is independent of ^23- Therefore, the S matrix elements 
obey relationships j23| 

Ser = Se^{c23 —S23, S23 C23), 
Stt = Sfj_f^{c23 —S23, S23 C23), 

Srfi = ~Sfj,T{c23 ^ —S23, S23 ^ 023). (14) 



B. Relations between neutrino oscillation amplitudes with NSI 



The secret behind the relations between Sei_i and Ser and the others in f|T^ is that H 
is independent of 6'23, or in other words, the invariance of H under the transformation 
C23 — > —S23, S23 C23. When NSI is introduced there exists the following additional term in 
H: 



H 



23 



^ee ^e/x ^er 

F* F F 

^e/i ^/^A* 

F* F* F 

^er ^ iiT ^rr J 



U23 = 



ee ^e/x ^er 

e/i ^/J-T 
F* F 



<^ee C23^efi — S23S er 

C23£l^ - S23£*eT C23£/xAt + Sss^rr " C23S23 (^^/xr + £^^r) ^23^ 



S23£efi + C23e er 



^23'^eii ~ ■'23Cer "-23'^MM ''23=^-"- " '-23''23V^Mr "r o^^j <-23^Mt ~ ^23^% + C23S23(^wt " ^ttJ 
•823^6// + C23^eT (^23^1lt ~ ^ 23^ iJ.T + C23S23(^/iAi ~ ^tt) S23^ + C^'^^^rr + C23"S23(^ait + ^^r) - 



(15) 
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Because H^^^ in (fT5l) does depend on 623, the S matrix relations as given in (fT^ do not 
hold. However, if we consider the extended transformation 



C23 ~^ ~S23, S23 C23, 



£^j,T ~^ £^ri ^fiT ~^ (16) 

it is easy to show that H^^^ is invariant under the transformation (fT6l) . It means that the 
5* matrix relations (IT^ hold even with NSI provided that we extend the transformation 
to the ones in fll6l) . It not only implies the existence of useful relations between the S 
matrix elements, but also serves as a powerful tool for consistency check of perturbative 
computation. We will see in Appendices. \M and [B] that the computed results of the 5* 
matrix elements, and hence the oscillation probabilities, do satisfy f|T^ with the extended 
transformation f|T6l) . 

As we will see in Sec. IVI Bt the invariance under the extended transformation f|T6l) entails 
a remarkable feature that the terms which depend on e's in the v^ — v-r sector in the oscillation 
probabilities P{v^ v^), P{i^^ —>■ i^t), and P(z/r ^t) are all equal up to sign. 



C. Phase reduction theorem 



Now, we present a general theorem on reduction of number of CP violating phases in 
system with NSI, which we call "phase reduction theorem" for short. By looking into 
the results of perturbative computation 28| it was observed that when the solar Am^i is 



switched off the oscillation probabilities with NSI depends on phases which come from NSI 
elements and 5 in a particular manner, e.g., |e|e*^'^"'"'^''. It was conjectured on physics ground 



that the property must hold in the exact expressions of the oscillation probabilities [31 



With vanishing Am^i the system becomes effectively two flavor and hence the observable 
CP violating phase must be unique. 

Here, we give a general proof of this property which is, in fact, very easy to do. We first 
notice a simple relation which holds in the absence of Am^i, 
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where A 



2E ' 



X 



6 + 



and UJ = S + 



[diag(e , 1, 1)]q,/3Z^/3, the evolution equation reads 



Then, if we use a new basis Ua 



dx 











= H 











(18) 
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It is obvious from (fT8|) that the system depends on only three phases x = S + (pefi, ^ = S + (f>eT 
and 4>fj,r out of four. This particular combination of phases is, of course, depends upon the 
specific parametrization of the MNS matrix. The phase factor attached to the transformation 
matrix in fll7l) does not affect the oscillation probability because it is an over-all phase. 

The similar treatment with the same transformation as in fll7p (diag.(l, l,e~*'^)) can be 
used to prove that the phase reduction occurs if 612 = (^23 = 0) even though Amg^ 7^ 0.^ 
If = it is obvious that there are no effect of 6. 

This completes a general proof that number of CP violating phases is reduced by one 
when the solar Amgi is switched off, or one of the mixing angles vanishes. We emphasize 
that this property has implications to the real world; For example, the phenomenon of phase 
reduction occurs at the magic baseline, |^ = vr, in the perturbative formula to be obtained 
in Sec. I VI I even though Am^i 7^ 0. 



V. PERTURBATION THEORY OF NEUTRINO OSCILLATION 
A. e Perturbation theory 

To formulate perturbation theory one has to specify the expansion parameters. We take 
the following dimensionless parameters as small expansion parameters and assume that they 
are of the same order:^ 

-T— ~ ~ ~ e {a,p = e,fx,T). (19) 

Whereas, we treat and ^^^^ as of order unity. We collectively denote order of magni- 
tude of the expansion parameters as e, and hence we call the perturbative framework the e 
perturbation theory. In the absence of NSI our formulas of oscillation probabilities, of course, 
reduces to the Cervera et al. formula 23] , which we call the SI second-order formula in this 



paper. Correspondingly, we call our second-order probability formula the "NSI second-order 
formula" . It appears that in the standard case this perturbative framework accommodates 
the situation of relatively large 61^ within the Chooz bound [i^, and applicable to wide 
variety of experimental settings. 

Another approach would be to just expand in terms of Sap which is assumed to be small 
without any correlation with other SI mixing parameters. If NSI elements are extremely 
small, much smaller than the SI expansion parameters, such first-order formulas of NSI 
would be sufficient. It would be the case of NSI search in the next generation experiments 
as discussed e.g., in 30l . 



On the contrary, it often occurs in deriving constraints on various NSI parameters that 
the bounds on the diagonal e's. See-, and Err, are sometimes milder than the ones on the 
off-diagonal e's by an order of magnitude. If it is the case, we may need to keep the higher 
order of the diagonal e's in e-perturbation theory, to e.g. e^, in probabilities to analyze such 
situations. We try not to enter into this problem in our present treatment. 



^ We thank Hiroshi Nunokawa for calling our attention to this feature. 

^ We do not take — S23 as an expansion parameter because a rather large range is currently allowed 
and the situation will not be changed even with the next generation experiments 24| . 
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B. Formulating perturbative framework 



We follow the standard perturbative formulation to calculate the 5* matrix and the neu- 
trino oscillation probabilities [3]. Yet, we present a simplified treatment which is suitable 
for higher order calculations. For convenience, we start by treating the system without NSI 
in this section. We use the tilde-basis z> = f/^ 
tilde-Hamiltonian is decomposed as H = Hq + Hi, where 



23 z/ with Hamiltonian H defined in flTT]) . The 



A 



rAix) 

1 



(20) 



Hi = A 



'13 





L cissise 



ci3Si3e-^^ 

-si. 



+ Ari 



where A = ^I^^^ ta = i^a{x) = Though our treatment can be easily generalized 

to cases with matter density variation as far as the adiabatic approximation holds, we 
present, for ease of presentation, the formulas with constant matter density approximation. 
To calculate S{L) we define VL{x) as 



„2 2 
''12'-13 

C12S12C13 
2 

- ~'5l2Cl3Si3e 



C12S12C13 
'^12 

-Ci2Si2Si3e 



-S?2Cl3Sl3e 
-Ci2Si2Si3e~'^ 



S-inS 



12"'13 



(21) 



iHox 



S{x). 



VL{x) obeys the evolution equation 



i-^fi(x) = HMx) 
ax 



where 



Hi = e^^«^^ie-*^o^ 
Then, Q{x) can be computed perturbatively as 



n{x) = 1 + (-i) / dx'Hiix') + (-i)2 / dx'Hiix') / dx"Hi{x") + 0(e^ 



(22) 



(23) 



(24) 



(25) 



where the "space-ordered" form in fl2^ is essential because of the highly nontrivial spatial 
dependence in Hi. Collecting the formulas the 5* matrix can be written as 



S{L) 



f/23e^'^»^r](L)[/i3 



Therefore, essentially we are left with perturbative computation of VL{x) with use of 
calculate the S matrix.^ 



(26) 
to 



^ Since Hi in (|2ip contains order terms in addition to order e terms the formal expression in (|25p includes 
terms higher than ©(e^) which are meant to be ignored. The same statement applies to the computation 
to be carried out in Sec. |Vll 
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C. Matter hesitation and unitarity 

One of the usefulness of the e perturbation theory is that it allows to prove the property 
"matter hesitation" , which is a characteristic feature of neutrino oscillation in matter without 
NSI in small ^13 regime. The matter hesitation refers to the property that the matter effect 
dependent terms in the neutrino oscillation probabilities P(z/q, — > z/^) = e,/i, r) are 

absent to first order in e. Namely, it hesitates to come in before computation goes to second 
order in e. Though its validity heavily relies on the particular perturbative framework we 
work in this paper, it explains why it is so difficult to detect the matter effect in many 
accelerator experiments. 

In fact, it is easy to observe the property of matter hesitation; It directly follows from the 
structure of the S matrix in (l26l) itself. _We first note that in the tilde-basis Hi is free from 
the mater effect and it exists only in Hq. Therefore, the matter effect dependence exists 
only in efi and er (and their conjugate) elements in Hi in f l24p . and they are of order e. The 
same statement follows for fl in (125|) . Then, the matter effect dependence in S" = e~*'^°^il 
is only in e/i, er and ee elements at first order in e. Notice that the final rotation in 23 
space to obtain the S matrix in (126!) does not alter the property of the S matrix. Therefore, 
no matter effect dependence appears in the oscillation probabilities to first order in e in ee, 
e/i — er and in the /ir sector channels for different reasons: In P(t'^ v^) and Pii"^ Vt) 
the matter effect is trivially absent to order e because of no dependence in the S matrix. 
In P{ve i^^t) and P{i'e ~^ ^t) it comes in only at order because the S matrix elements 
are of order e. In Piye ^e) the matter effect is absent to order e because it is contained 
in a phase factor of the S matrix element. This completes the derivation of the matter 
hesitation, the property that matter effects comes in into the oscillation probability only at 
second order in e. 

We stress that the absence of the matter effect in the oscillation probability to first order 
in e is highly nontrivial, in particular in See- Since the matter effect coefficient va = 
is zeroth order in e it can affect the 5* matrix in all orders of e. But, in fact, absence of 
matter effect to first order in e in P(t'e ^e) can be understood by unitarity. For the most 
nontrivial channel, the relevant unitarity relation is 

1 - P{Ve ^ Ve) = P{Ve ^ 1^^.) + P{l^e ^ ^r) (27) 

Since P{ve ^fi) and P{ve J^r) are at least of order as will be shown in Sec. I VI I and 
Appendix [HI the matter dependent term in P{ve ^e)-, which is involved in the left-hand- 
side in (1271) . has to be second order, or higher, in e. It should be noticed that this argument 
is valid not only in systems with SI only but also in the one with NSI, the matter hesitation 
property for P{ve —>■ ^e) in the presence of NSI. Also notice that the same argument does 
not go through for 1 — P(z/^ — > v^) because P^i^^ Vt) can contain the terms lower than 
e^. If fact, there exists the first order term in e in P{v^ v^) and -P(z^^ Vt) which are 
proportional to the matter effect coefficient a. 

D. Implication of matter hesitation to neutrino oscillation with NSI 

There is a clear implication of the property of matter hesitation to the system with NSI; 
The terms with the NSI element See must appear in the oscillation probability only at third- 
order in e or higher. It is due to the special nature of Eee that can be introduced as a 
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renormalization factor of the matter effect coefficient a, a — a(l + Eee)- Since the terms 
with a are aheady of order e^, the terms with See must be at least of order e^. 

The reader should be puzzled by the above statement. One may argue quite naturally 
that there must exist a term with first order in See in the survival probability P(fe ~^ ^'e)- 
In fact, such a term does exist in the relevant S matrix element as one can see in (lAip : 

See = e-'''^^^|l - tisl^TA + £eerA)AL} (28) 

Resolution of the puzzle, therefore, is that the first order term of See cannot appears in the 
oscillation probability because it is purely imaginary, or a phase ignoring terms. However, 
to confirm the cancellation of second order term we must go beyond the present treatment 
by keeping the order terms in the S matrix. It will be done in the next section. 

A more general question is whether the matter hesitation can be generalized into the 
whole systems with NSI. We have already answered the question at the end of the previous 
subsection. This feature is to be verified by explicit computation in Sec. |Vll 



VI. NSI SECOND-ORDER PROBABILITY FORMULAS 

Now, we present the expressions of the oscillation probabilities with NSI which is valid 
to second order in e. For ease of computation we use a slightly different basis which we call 
the double-tilde basis with Hamiltonian 



(29) 



and the corresponding S matrix 



s{L) = U2zUiAl)uIuI, 



23 



(30) 



where S{L) = Texp —i dxH{x) . The zeroth order and the perturbed part of the 
reduced Hamiltonian H are given by 



A 





1 



^12 





C12S12 

C?2 





"^13 


L Ci3Si3e' 



ci3Si3e-^^ 

sj. 



+ AtaU 











e* 






f/l3 


e* 


e* 


^TT _ 





(31) 



(32) 



where A 



^"^31' 



2^ , ta = ^^i"- , ta = ■ To simplify the expressions of the S matrix 

elements we use the NSI elements in the tilde basis, Eap = {U^^) a'y£'ys{U23) , defined in ([T^ . 
Notice that e's are invariant under the extended transformation flTBl) . 
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The perturbative computation of the S matrix elements can be done with the formulas 
similar to the ones in the tilde basis in Sec. IV B[ In this section we concentrate on the struc- 
tural analysis of the NSI second-order oscillation probabilities e.g., for analysis of parameter 
determination. We collect all the resultant explicit expressions of S matrix elements and 
the oscillation probabilities in Appendix [X] and [B], respectively. The results of third-order 
calculation which are necessary to complete Table [I] are presented in Appendix O 



A. Electron neutrino sector 

The most distinctive feature of the NSI second-order oscillation probabilities in the z/g- 
related sector is that they have very similar forms as the SI second-order formulas 2^ but 
with the generalized atmospheric and the solar variables: 

a a 

E = ( Ci2Si2^^ + C23£e/. " S23£er ) = Ci2Sl2^^e'^ + \€e^,\e'^^^ (33) 

\ Qj J d 

and their antineutrino versions 

- _ ^"^31 I M r pMp-*'^ - ICQ ^"^31 I l~ \^-i4>er 

^± = —SrS r [S23£f^ij_ -\- C23£f,T-)e — +S13 h l^erie , 

a ^ a 

-ci2Si2^^ + 023^:^ - S23e:r) e"^' = -012^12^^6-^^ + \ee^\e-'^^- , (34) 

where c^ea = i^ + ^ea (tt = yU, t). The particular dependence on NSI elements in fl33l) and flM|) 
has root in the form of the perturbed Hamiltonian f[T^ in the tilde-basis, from which it can 
be understood that e^e/x ^ind e^t play the role of the mixing angles which govern 1-2 and 1-3 
transitions, respectively. At the second equality in the right-hand-side of these equations we 
have introduced a new notation Amg^^ = -^bm^^ where ± sign indicates the sign of Amg^^, the 
mass hierarchy, and 5m|j^ = |Am|^| > 0. Note that a = 2\/2GFNf.E > 0. For convenient 
notation we parametrize these quantities as 

e± = |e±|e*^±, S = |S|e'«. 

e± = |e±|e*^'±, S = \E\e'^. (35) 
To represent the oscillation probability in a compact way we define 



o sm^ — — . 



Y+ = sm — sm — — L 



^ \5ml^TaJ AE"'" AE 

Z^sin^^. (36) 



AE 

For anti-neutrinos we have flipped sign of a, and hence 



Y± = Y^. (37) 
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Z is obviously invariant under the sign change of a, Z = 

With these notations and by defining A31 = for 
oscillation probabilities P(z/e Ue), P{^e i^fi) and P(i^e - 
neutrino counterparts of the latter two) can be written as 



Pi^, 
P{l^e 
P{l^e 

Piv., 



1 -4x±|e±|2 



44^±|e±p + 44z|sp + 
44x±|e±p + 44z|s|2- 

44^±|e±p + 44z|sp + 

CP[P(z/e^z/^)] 

44^Ti0±i'+44^isr+ 

T[P(Z7, ^ Z7^)] = TCP[P(i/e 

44^t|0±I' + 44Z|S|2 + 



BC23S23>±| 
3C23S23Kt| 

BC23S23^±| 
BC23S23Kp| 

3C23S23Kpl 
BC23S23Kf| 



simplicity of expressions, the 
-> Vt) (together with the anti- 



e± 

6+ 


cos(^ - 
cos(^ — 


6+- - 


IA31I) 

1 ^.51 1 / 


(38) 
(39) 
(40) 


e± 


cos(^ - 


e± + 


IA31I) 


(41) 


e± 


cos(^ - 


e±- 


IA31I) 


(42) 


e± 


cos(^ — 


e±- 


IA31I) 


(43) 


e± 


cos(^ - 


e± + 


IA31I) 


(44) 



The upper and the lower signs in the above equations are for the normal and the inverted 
hierarchies, respectively. The expression of P(z/e — > i^e) is so simple because of the unitarity. 

Notice that to second order in e, the oscillation probabilities in the related sector 
do not contain any NSI elements in the z/^ — z/^ sector, e^r etc. It should not come as a 
surprise because in the z/^ and appearance channel from the leading term of the 5 
matrix is already of order e, and it can contain only the Ue related NSI elements, Sen and 
Est- Therefore, to order there is no room for NSI elements in the z/^ — Ur sector in the 
appearance probabilities. We will see in the next subsection that this simple fact leads to a 
great simplification of the oscillation probabilities in the z/^ — Ur sector. 

We note, in passing, that because of the relation Y± = \/X±Z which is easily recognized 
by dSS]) it is evident that the oscillation probabilities can be written in a form of absolute 



square of addition of the atmospheric and the solar terms. ^ For example, Pive 
the form 



takes 



PiVp 



+ C23V^|5|exp \i {i-Q±- IA31I 



(45) 



At the magic baseline, = vr, 
simple expression of the oscillation probability, P{y^ 



the second term vanishes because Z 

v^) =4s23X±|e±| 



0, leaving a very 



B. Vyi — Vt sector 

As in the z/g-related sector there is a distinct characteristic feature of the oscillation 
probabilities with NSI in the — sector. Namely, to second order in e they can be 



We thank Stephen Parke for caUing our attention to this point. 
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decomposed into the three pieces with different dependences on NSI elements, the vacuum 
term, the ones with Sais in the i^e-related and the z/^i — z/^ sectors, respectively: 

P(z/„ Uj3; Eeij,, Ber, £(1(1, £,it, £tt) = -P(«^q i^iB] 2 ffavor in vacuum) 

+ P{i^a ^ i^p]£f,f„£f,T,£TT) (46) 

where a and (3 denote one of fi and r. The explicit expressions of these terms will be 
displayed in Appendix [Bl 

The point is that the last term in fHBl) is universal, up to sign, among all the three 
channels, Pii^^ Vfj)-, P{^fi ~^ ^r), and P(z/t- —>■ 1/^). Though it may look mysterious, it is 
in fact very simple to understand it. By unitarity it follows that 

P{^tM ^t,) + P{^t, ^Vr) = 1 - P{i^^ z/e), 

P{Ur ^ Ur) + P{Ur ^ U^) = 1 - P(z/, ^ Z/^). (47) 

We note that P{i^^ — ^ z/g) and P^Ur ^e) do not contain e^^. Err, and e^,- to second 
order in e. Then, it follows from the first equation in fHTj) that Pi^UpL — > i^T',£fi^i,£iiT,£TT) = 
—P{Ufj_ — > z/^; e^/i, e^xr, ^rr)- Notlclug that the terms related to e's in the z/^ — Ur sector are 
T-invariant, the relations P(z/t- J-'t]£^i^ii£^iti£tt) = ~P{^fi ~^ ^T]£^J.^i,£^lT,£TT) must also 
hold. Therefore, the e^^g (a, /3 = /i, r) dependent term in the three channels are all equal up 
to sign. The equality P{iy^ t'/.; ^m^, e^/^r, e^rr) = P{i^t ^ l^T^,£^,^^,£^,T,£TT) also follows from 
the relationship between S matrix elements due to the extended transformation fll6p . 

VII. PARAMETER DETERMINATION IN NEUTRINO OSCILLATION WITH 
NSI 

Thanks to the NSI second-order probability formulas derived in the previous section, 
we can now address the question of how simultaneous measurement of the SI and the NSI 
parameters can be carried out. However, we must first warn the readers that our discussions 
in this section are based solely on the NSI second-order formulas, and hence its validity 
may be limited. Nonetheless, we believe that ignoring the effects is quite safe because we 
anticipate e ~ 10~^ in our perturbative framework. 

A. SI-NSI confusion 

One of the most distinctive features of the oscillation probability formulas in Sec. |Vl] is 
that the NSI parameters e^a = /i, t) appears in the particular combination with the SI 
parameters as in (1551) and ([31]). What that means in the context of parameter determination? 
It means that, in general, determination of SI mixing parameters, ^13 and 6, has severe 
confusion with determination of NSI parameters Eap, and vice versa. However, it should be 
noticed that it does not mean something like "No Go" theorem. Namely, there is a way to 
circumvent this problem. It is a complete determination of the SI and the NSI parameters, 
the possibility we address later in this section. 

Nonetheless, we should note the following: If such complete determination is somehow 
not feasible experimentally, our result may be interpreted as an analytic proof of the "NSI-SI 
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confusion theorem".^ It is a powerful statement because it not only reveals the existence of 
confusion but also illuminates which SI parameters are confused with which NSI parameters 
via which manner. 

In fact, the characteristic feature in (I33p . namely, ^e/x only couples to the solar scale 
oscillation and ier the atmospheric one, would affect the resolution of the ^13-NSI and the 
two-phase confusions. Coupling between the solar and the atmospheric degrees of freedom 
bridged by a NSI element is the key to the resolution of the 6'13-NSI confusion by the 
two-detector method [3l[. Therefore, the resolution mechanism might be affected by the 
simultaneous presence of two e's, which "decouples" the solar and the atmospheric degrees 
of freedom. This point deserves a careful investigation. 



B. Strategy for parameter determination 

To gain a hint of how we can proceed let us look at Table. [H We first note that it is not 
possible to detect the effects of See because it is of third order in all channels, and hence we 
have to omit it from our subsequent discussions.^'^ It is also well known and is obvious from 
the probability formulas in Appendix [B] that and Srr come in through the form — Err 
and therefore only their difference is measurable. 

Next, we observe that in i/^ and z/^ appearance channels from Ug, only the u^, related NSI, 
Eefj, and Eer appear to second order in e. Therefore, the obvious strategy is to use these 
channels for complete determination of them simultaneously with and 6. Then, we may 
be able to determine the rest of the NSI parameters in the — z/^ sector by disappearance 
and appearance measurement in that sector. 

The important point is therefore that one can explore the effects of Ee^ and Eer in Ue related 
channels while ignoring e^^, and Err- It is a good news because the appearance channels, 
assuming excellent detection capability of z/^ and u^-, have great potential of detecting the 
effects of NSI [311 • Once Ee^ and Eer are measured one can proceed to determine the rest of 



the NSI elements and e^t using the oscillation probabilities in the — Vr sector. 



C. Complete measurement of the SI and the NSI parameters; ^13, 6, and Eer 

Now, we start to formulate a recipe for complete determination of the SI and the NSI 
parameters. Based on consideration in the previous subsection, we concentrate on P{v^ 
u^) and P(z/e — ^ z^r) and their CP and T conjugates. By looking into the expressions of 
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We note that a different type of the confusion theorem was derived in 27[ which involves 6*13 and NSI 



parameters in production and in propagation processes that obey a special relationship. 
If we take the setting with only See as NSI, it can be regarded as uncertainty in the matter density and 



it is known that neutrino factory has a great sensitivity to it |47l . I48j . However, in our current setting 
the issue of matter density uncertainty is much more severe and universal; It produces uncertainties in 
determining all the NSI elements. Clearly, the discussion of this point is beyond the scope of the present 
paper. 

If 6*23 is deviated significantly from the maximal so that cos 26*23 ^ ^, then the terms with e^^ — Ett can 
have sizes of order e. In this case, it may be possible to detect the effects of e^^ — Err and measure (or 
constrain) it even without having a priori knowledges of Eefi and Eer- 
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oscillation probabilities in (139!) and (HOl) (and other related ones which will be given below) 
one notices that the observable quantities are of the forms 

|9±|^, — 9±, in neutrino sector, and 

|6-i-P, 9±, in antineutrino sector. (48) 

where the phase ^, 9±, etc are defined in There are altogether six quantities. 

Suppose now that somehow we were able to determine all these quantities. We discuss 
in the following subsections how it can be done. Here, we show how they determine the SI 
and the NSI parameters, S13, 6, \ee^\, je'erl, 4>efi, and (per- It niay be sufficient, assuming that 
the inversion is possible, to express the observable in terms of the physical parameters. We 
start with the neutrino sector: 



a 

IC12S12 +\^eij,\ +'2^Ci2Si2\eeA COs(() - 0e/.) • (49) 



For phase difference we obtain 



Pi* " 1 



C12SV 



X 



By taking the real and the imaginary parts of flSUj) one can obtain cos 2(^ — 9±) and sin 2(^ — 
9±), respectively. For antineutrinos we obtain 

+ |^erPT2Si3|eer| ( j COS0er, 

- 2 / ^"^21 \ ^ i~ i2 i'^"^21 /r " \ 

FI = I C12S12 I + l^e/xl - 2Ci2Si2|ee;,| COS(()-0eM)- (51) 

Similarly, the equation for the phase difference 9± — similar to (I3U]) can be obtained by 
making the transformation a — > —a, 6 — > —6, (pe^ — > — 0e^i, 0er — > — 0er in fl50|) . 

Having the six equations altogether with given six observable, |6±|, |0±|, ^~9±, 
and ^ — 9±, they can be solved for the six unknowns, S13, S, two complex numbers iefj., and 
^er- Given the latter two numbers one can determine the original and Eer- Therefore, 
the rest of the problem in simultaneous determination of the SI and the NSI parameters is 
how to measure the above six observable. 
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D. Measurement with a monochromatic neutrino beam; fe sector 



In this subsection, we discuss a way of determining the SI-NSI combined parameters in 
l8|) by assuming a set of measurement at an energy E, aiming at their complete determi- 
nation.^^ Though it might not be a practical way, by describing a concrete method we try 
to illuminate characteristic features of the problem of complete determination. With the six 
unknowns we have to prepare neutrino oscillation measurement of six different channels. 
Suppose that one measures the following six probabilities at a neutrino energy E, P{ve — >■ 

z/t-) = CP[P(z/e ^t)]-, P{i^iJ. ^e) = T[P(z/e z/^)]. Noticc that we have intensionally 
avoided to use the channels which require Ur beam which, if not impossible, would be very 
difficult to prepare. From (139|) . (140!) . and fHTl) . it is easy to obtain 

Piiy,^iy^) + P{iy^^iy,) = 8sl,X±\Q±\^ + 8cl^Z\E\^ 

+ 16c23S23l±|H||e±|cos(e-e±)cos|A3i| (52) 
P(z/e z/^) - P(z/^ z/e) = 16c23S23Kt | H 1 1 Gi | sin(^ - 6^^) siu | A31 1 (53) 

Pil^e + P{ye Vr) = 4X±|e±P + AZ\E\^ (54) 

Similarly, for the antineutrino channels, we obtain from fH21) . fl43p . and 



P{ue^l^^.) + P{l^^.^l^e) = 84X^|e±|2 + 84Z|S|2 

+ 16C23S23'i^Tl'='l|0±l COs(^- ^±) cos IA31I (55) 

P(z?e z/^) - P(z/^ z/g) = 16c23S23KfI^I|0±I sm(^- ^±) slu IA31I (56) 

P{l^e ^ + ^('^e ^ ^^r) = ^X^\Q±\^ + ^Z\E\^ (57) 

It is easy to solve these equations to obtain |6±|, and (^ — 6'-|-) (for neutrinos), and |0±|, 
and — 9±) (for antineutrinos). 
It may be obvious that the above analysis can be converted to the rate only analysis by 
replacing the probabilities P(t'Q, — > z/^) by energy integrated number of events with fluxes 
and cross sections / dEF^ayj'qPii'a ^fs), and the similar integrated quantities of X± etc. 



E. Determining the NSI parameters in the v^j^ — sector 

After measurement of 6*13, 5, and and Ser as described in the previous subsection, 
one can proceed to determination of the NSI parameters in the — sector with (for 
concreteness) mono-energetic beam. As we saw in Sec. I VI BI the and e^r dependent term 
in the oscillation probabilities is universal in P(z^^ — > z/^), P{y^ —>■ i^t), and P(z/t- — > Ur). 
Therefore, one can simply use one of the above three channels, which means that r neutrino 
beam, even if it were prepared, does not help. 

The oscillation probability P(z^^ z/^) derived in Sec. |Vl]can be written as 

P{y^i ^fjLi ^efj,: ^er; ^fi^iy ^/^r; ^tt) P^^fjb ^ ^efiy ^er) 

+ S^'\e,, - Err? + W(°)(^MM - err)Re{e^,) + Q(°)Re(e^.)2 + J(°)lm(e^,)2 (58) 
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It was proposed that such a monochromatic neutrino beam can be prepared for and Ve beams 49|, |50| • 
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where the exphcit form of the coefficients can be easily read off from the expressions in 
Appendix [B] and we have the similar expression for antineutrinos. We have obtained two 
equations for the three unknowns, — £rT, ^e{e^r), and Im(e^^). Clearly we need one more 
equation to determine the three unknowns, which is unavailable under the current setting. 
Thus, we have to conclude that a complete determination of the NSI elements in the — Uj- 
sector is not possible by measurement at a monochromatic beam or the rate only analysis. 



F. Necessity of spectrum analysis 

Doing measurement at six different channels is not the unique way of carrying out com- 
plete determination of the six parameters. Even in the case where only the "golden channel" , 
P(z/e Vfj) and P{i'e ^'/i), is available, one can in principle determine |0±P, ^ — 9±, 
and their antineutrino counterparts by spectrum analysis. It is because the energy and base- 
line dependences of the coefficients of these quantities in the oscillation probabilities in fl39|) 
and (j^2l) are different with each other. In the — sector all the NSI elements cannot be 
determined by the rate only analysis, and need for the spectrum information is mandatory 
in this sector. 

It appears that one of the most promising ways to carry this out is the two-detector 
method [sij. It has been applied to the Tokai-to-Kamioka-Korea (T2KK) two-detector 
complex which receives neutrino beam from J-PARC 42, 52, 5^|.^^ In the context of neu- 



trino parameter determination in neutrino factory with NSI as well as SI, this method was 
examined in detail in 31 



G. Parameter degeneracy; Old and new 

1. NSI- enriched conventional type degeneracy 

The parameter degeneracy is the problem of multiple solutions in determination of lepton 
mixing parameters |33l. l34j . l35l |. It is known to be a notorious problem for their precision 
measurement. See (sd. IstI] for a global overview of the degeneracy, and (stI . 58| for pictorial 



representation. 

We give evidences that the phenomenon has an extension to the system with NSI.^^ Our 
discussion based on the matter perturbation theory in Sec. I Villi indicates that the parameter 
degeneracy prevails in system with NSI but with new form which involve NSI parameters. 
Set of equations for observable we have derived in Sees. IVII CI and IVIIDI shows that the 
sign-Am^i and the 6*23 octant degeneracies exist because the equations take different form 
for different mass hierarchies and octant for a given set of observable. It is also very likely 



Other possibility would be the one called the "on axis wide-band beam approach" which was proposed in 



a concrete form in the project description for Brookhaven National Laboratory [54|. Precise estimation 
of the potential in doing spectrum analysis, however, depends upon which kind of detector is chosen and 
the ailtiial performance of the detector. 

See 551 for effects of the systematic errors and optimization of the similar two-detector setting in param- 
eter determination in neutrino factory. 
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Notice that introduction of NSI parameters leads to a new solution of the solar neutrino problem [59 1 . 
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that the intrinsic-type degeneracy survives with a NSI-enriched form, as one can see in the 



bi-probabihty diagram [34| given in Fig. 2 of [31 



2. New type of degeneracy 

Here, we present a completely new type of parameter degeneracy which may be called as 
the "atmospheric-solar variable exchange" degeneracy. We work with the setting of measure- 
ment of six channels at a monochromatic energy. First of all, one notices that determination 
of the neutrino (un-barred) and the antineutrino (over-barred) variables decouples with 
each other. We discuss only the neutrino variables below because the antineutrino ones is 
so similar. To simplify the expressions we restrict ourselves to the case of maximal 623- By 
combining (IH^ . (IH^ . and it is easy to show that the phase variable can be determined 
as 



tan(^ - e±) = cot |A 
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P - P 



(59) 



where we have used a simplified notation P^ji = Piya ^13) ■ 

It is easy to show that if the mass hierarchy is known the solution of this equation is 
unique in the physical region — vr < ^ < vr and —it < 6± < tt. Then, the solutions for |B±| 
and |S| are given by 



p + p 

P + P 
8Z 



1± Wl- 



sin^(e-e±)sin^A3i \Pe^ + P, 



P — P 



ITA 1 



sin^(e-^±)sin^A3i \Pe^. + P. 



P — P 



(60) 



It is the solar-atmospheric variable exchange 
G^'*! and IS'^^-'I, then the second solution 



Notice that the degeneracy is quite new; 
degeneracy. That is, if there is a solutiun 

IG^"*! = and = \J ^|g£'*| exists. Notice that the new degeneracy does 

not survive when NSI is switched off where ^ = 5 and 9± = {). Namely, there is no phase de- 
gree of freedom in the atmospheric variable in the limit, while only phase degree of freedom 
exists in the solar variable. 

Now we turn to the sign-Am^ degeneracy. At first sight there is no sign-Am^ degeneracy 
because the sign-Am^ flipped solution of ^ — ^± has to satisfy the same equation (!59|) which 
has no explicit dependence on the sign. Nevertheless, there is indeed a sign-Am^ flipped 
solution. If ^^^^ and are the solution to (159!) then there are another solutions ^^^^ = ^^^-'ivr 

and 9^^^ = 6^}^ =p TT. It means the existence of the sign- flipped solution of G and S, which 
can be another solutions if accompanied by (Amg]^)*^^-* = — (Amg^)^^^ With these solutions 
of the phase equation there exist the similar degenerate solutions as in fl60|) . Again the 
sign-Am^ degeneracy does not survive in the no NSI limit because of no degrees of freedom 
of 6± in the limit. In conclusion we have uncovered new degeneracies of the intrinsic and 
the sign-Am^ flipped type which exist as a consequence of the presence of NSI. 
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VIII. MATTER PERTURBATION THEORY WITH NSI 



As a first step toward understanding the degeneracy we examine neutrino oscillation with 

a. 



It is known 34 that 



NSI by matter perturbation theory following the treatment in 
structure of parameter degeneracy is particularly transparent in the region where the matter 
effect can be treated as a perturbation, as explicitly verified in the analyses in 42. 5^ 
for further explanation of this point. 



See 



For simplicity, we restrict our discussion to Vf. related appearance measurement in this 
section. In concordance to these works we consider v^. and z/g appearance measurement with 
conventional muon neutrino beam and its antiparticles. 



A. Structure of the oscillation probability with NSI in matter perturbation theory 

If we restrict ourselves into the first order in a, the matter effect coefficient, the only 
terms that survive are the ones up to first order in ^e/x oi' ^er- The oscillation probability in 
channel is given to first order in matter perturbation theory as 



(61) 



where the leading term is the Arafune-Koike-Sato (AKS) formula without NSI [18 
P{u^ ^u,;e = 0)aks = sin^ 2^135^3 sin^ A31 + cj^ sin^ 2^12 y-^^- J ^ 
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4l 



+ 4X 



A31 [cos 6 sin 2 A31 — 2 sin 6 sin^ A31] 



'23 \^4^ 



+ 2 sin^ 2^133^3 - 



1 2 A 1 
— sm A31 - - sm 2 A31 

^31 ^ 



(62) 



In 



, ^31 



4E 



I = 2\/2GFNeE as before. Jr (= 012^120^3513023523) denotes the 
reduced Jarlskog factor. 

The first order matter corrections which include the first order NSI effects in e's can be 
obtained by taking the first order term in a as 



P(iy„ 



T-'e] ^erjNSI 



ae) 



X 



C23S23'5l3 \ kerl COs(5 + 0e 



sin^ A31 1 



31 



sin 2A31 1 + l^er I sin((5 + (per) sin^ A 



31 



- C12S12C23S23 



Am2^ 



l^erl cosher I A31 - - siu 2 A31 ) - l^er | siu siu A31 



, (63) 
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We got rid of a higher order term which was kept in our previous references, e.g., 34 l42l |52|. 
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X 



^e] ^ea)NSI 



C0S(5 + 0e/x) 



aL \ 



'23 



A 



31 



- ^sin2A3i^ +C23|ee/.|sin(5 + 0eM)sin^^3i| 



- C12S12C23 



Am^i 



£efj.\ COS( 



C23A31 + — sin2A3i ) + S23\eef^\ sin sin A31 



(64) 



The antineutrino probability P(z^^ Ve'-, ^eri ^ep) can be obtained by making the replacement 

Notice that both of the CP violating leptonic 



-0/3- 



in fl6T|) : a —a, 5 — > 27r — 5. Ea/s 

KM phase 6 and (pais due to NSI elements changes sign when we discuss the time reversal 
process t'^ — > z/g, as opposed to z/g — > z/^ in the previous sections. 

For the oscillation probabilities in the z/^ — z^r sector we only deal with the one in z/^ 
disappearance channel (which may be easiest to measure) to first order in e: 



P{^ii ^fj.] 1st order in e) 



1 - 4C23S23 sin^ A31 + 4c?,4^s? 



'12'-23''23 



Am|i 



+ 2C23S23 



('^23 ~ -^23) (^MM ~ ^tt) — 4c23S23Re(£:^T- 



A31 sin2A3i 
aL 



2E 



sin 2 A31 



— 8C23S23(C23 — S23) 



C23'523(^u/t ~ ^rr) + [C] 
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5^3) Re (e 



sin2A3i. (65) 



Notice that (l65l) is already in the form of first-order formula in matter perturbation theory. 



B. Sign-Am^ and ^23 octant degeneracies prevail in the presence of NSI 

In this subsection, we discuss the fate of the sign-Am^ and the ^23 octant degeneracies 
in the presence of NSI. In the conventional cases without NSI, they are known as notorious 
ones among the three types of degeneracies because they are hard to resolve and the former 
can confuse CP violation with CP conservation. The sign-Am^ degeneracy was uncovered 
in systems without NSI by noticing that the oscillation probability P(z/^ — > z/g) in vacuum is 
invariant under the transformation Am|^ — > — Am|]^, 6 ^ tt — 6 without changing ^13 34 



It maps a positive Am|^ solution to the negative one, and vice versa. The presence of the 
symmetry as well as the fact that it is broken by the first order matter terms can be seen in 



Now, we observe that the sign-Am degeneracy prevails in the presence of NSI. That is, 
the NSI induced terms in the probability (jHSD and flM|) . though they are "matter terms", 
are invariant under the extended transformation 

Amgi -Am^i, 

6 TT — 6, 

(pea ^ 27r - 0e„. (66) 

while keeping ^13 and \eea\ fixed, where a = /i,r.^'' The symmetry is broken only by the 



Under the transformation ((66)) . the trigonometric factors in ((64|) and ([63)l transform as foUows: cos{S + 

4>ea ) — C0S((5 + 4>ea ) , Sin{S + (pea ) + Sin((5 + (j)ea ) , COS (j)ea ~^ + COS (pea , and sin (j)ea — sin ( 
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matter term in (162|) which is independent of NSI; The symmetry is broken by the matter 
effect which has exactly the same magnitude in systems with and without NSI. Therefore, 
to first order in matter perturbation theory, the sign-Am^ degeneracy exists in systems with 
NSI to the same extent as it does in the system without NSI. Given the robustness of the 
sign-Am^ degeneracy in the conventional case we suspect that the degeneracy in systems 
with NSI has the similar robustness. 

Similarly, one can easily show that the ^23 octant degeneracy survives the presence of 
NSI. It can be readily observed that -P(z^^ — > t'^; 1st order in e) in (165]) is invariant under 
the transformation 



C23 S23, 
S23 ^ C23, 

{£,jifi - £rr) -{^fifi - ^tt)- (67) 

It means that the 6'23 octant degeneracy prevails in the presence of NSI, and actually in an 
extended form which involves NSI parameter — Err- Since this NSI parameter decouples 
from Pii^fj, Ve) to second-order in e, the presence of the 6'23 octant degeneracy remains 
intact when the NSI is included though values of the degenerate solutions themselves are 
affected by the presence of E^a- 

It is interesting to note that both of the two degeneracies discussed in this subsection 
have common features. Their presence can be discussed based on (approximate) invariance 
under some discrete transformations, and with NSI the transformations are extended to the 
ones which involve NSI parameters. Most probably, our treatment here is the first one to 
signal the existence of the degenerate solutions which involves both the SI (6'i3 and 5) and 
the NSI parameters. 



C. Decoupling between the degeneracies in the presence of NSI 



In [42] the property called "decoupling between degeneracies" are shown to exist for 
experimental settings with baseline shorter than ~ 1000 km which may allow treatment 
based on matter perturbation theory. See also 3^ and 61 1 for preliminary discussions. 



The property of decoupling between degeneracies A and B guarantees that when one tries 
to resolve the degeneracy A one can forget about the presence of the degeneracy B, and 
vice versa. Existence of NSI terms, in general, influences the discussion of decoupling. It 
is the purpose of this and the next subsections to fully discuss the fate of the decoupling 
in the presence of NSI. Since it is one of the most significant characteristic features of 
the degeneracies in matter perturbative regime, we believe it worth to present a complete 
treatment. 



1. Definition of decoupling between degeneracies 



To define the concept of decoupling between degeneracies A and B, we introduce, following 
4^ . the probability difference 



- P - z.,; (AmL)(^ e^l efl 4")) , (68) 
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where the superscripts a and b label the degenerate solutions. Suppose that we are discussing 
the degeneracy A. The decoupling between the degeneracies A and B holds if AP"* defined 
in (l68l) for the degeneracy A is invariant under the replacement of the mixing parameters 
corresponding to the degeneracy B, and vice versa. 



2. Matter-perturbative treatment of the degenerate solutions 

We follow 4^ to define the degenerate solutions in a perturbative manner. Throughout 
the discussion in this section we assume that deviation of 623 from the maximal angle 7r/4 
is small. A disappearance measurement, — *■ f^, determines S23 to first order in 5^3 as 
{^23^^^ = {-^23)^^^^ + -^is)' where {s^^)^^^ is the solution obtained by ignoring sf-^. It is given 
by (S23) = ^ 1 ± a/1 — sin^ 2^23 • In leading order the relationship between the first and 



the second octant solutions of 623 is given by s^f = (3^^. 

A Vf. appearance measurement determines the combination s^s sin^ 2^13. The first and the 
second octant solutions of Q23 are also related to leading order by ^23*5^3* = s^^s^^'^ . In an 
environment where the vacuum oscillation approximation applies the solutions corresponding 
to the intrinsic degeneracy are given in Appendix [D] as 




■ r '^13 • r 

smd2 = — ^smdi 

^13 



where 



cos h = (e^^ COS 5i + ^ ) (69) 

^13 ^ 



Y 

= sin 2^12 cot ^23 A21 cot A31 . (70) 

X 



and the superscripts (1) and (2) label the solutions due to the intrinsic degeneracy. The sign 

(2) (2) 
=1= for cos52 are for = ±|5^c|; and in the solution of 6 is meant to be the ^^3 solution 

given in the first line in (!69|) . 

As we saw in the previous section, an extended form of the sign-Am^ degeneracy is given 
under the same approximation (mod. 27r) as 

/inorm flinv f /\^2 N norm f \ mv rnorm xmv (i N norm (i \inv (^'\\ 



More precise meaning of the term "perturbative" is as follows: Since the disappearance probability by 
which 023 is determined is of order unity we disregard quantities of order e or higher. They include the 
matter effect, 6*13, and NSI. Similarly, I'e appearance probability is of order the relationship between 
the two degenerate solution inevitably contains a small quantity, which is 613 in this case. But, all the 
quantities of higher order are neglected. If the near-far two detectors are involved, like in the case of 



T2KK [4^|5^, the degenerate solutions are essentially defined by the near detector. In this case, the 
second detector is meant to give raise to perturbation effect to lift the degeneracy. For more concrete 



example of this feature, see [42 1 
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where the superscripts "norm" and "inv" label the solutions with the positive and the neg- 
ative sign of Am|]^, and (j)ap denotes the phase of Sap- The validity of these approximate 
relationshi ps i n the actual experimental setup in the T2K II measurement is explicitly veri- 



fied in j42|, |52| . It should be noticed that even if sizable matter effect is present the relation 
(I7T]) holds in a good approximation if the energy is tuned to the one corresponding to the 
vacuum oscillation maximum, or more precisely, the shrunk ellipse limit [g^ . 



D. Decoupling between the sign-Am^ and the ^23 octant degeneracies 



Let us start by treating the sign-Am^ degeneracy. For this purpose, we calculate 



norm mv/ 



as defined in 



Apnormmvj^^^ _^ ^^-j defined in (1681) . Thanks to the extended symmetry (IMj) obeyed 
by the appearance probability, it is given by the same result obtained without NSI in 42 |: 



AP 



norm mv/ 



sin^ 29 



norm / „ norm N 2 



13 



'23 



(A 



31 j 



■sin^(A, 



31 y 



sin2(A3i: 



(72) 



where the superscripts "norm" and "inv" can be exchanged if one want to start from the 
inverted hierarchy. Therefore, breaking the sign-Am^ degeneracy requires the matter effect 
but not more than that required in resolving it in systems without NSI; NSI does not 
contribute resolution of the sign-Am^ degeneracy but it does not add more difficulties. 

By following the same discussion as in [42] , we observe that AP ™^ is invariant under 
the transformation 9^1^ ^ 9^^'^ and 9^1^ ^ 9^^'^, because AP depends upon ^13 

and ^23 only through the combination sin^ 2^13533 within our approximation. Therefore, 
resolution of the sign-Am|^ can be done in the presence of the ^23 octant degeneracy. 

What is the influence of the i/^ disappearance channel in the discussion of decoupling? 



Using the first-order formula in (!65|) . AP 



norm mv/ 



can be computed as 



AP 



norm mv 



'^12^23 -523 



+ Ac: 



23 "^23 



--23 



'23 



Aml^L 



sm 



Aml^L 



4E 

4c23S23Re(e^^ 



2E 
aL 
2E 



sm 



Aml^L 
2E 



(73) 



It is manifestly invariant under that the transformation in fl67|) . and hence the sign-Am^i 
degeneracy decouples from the 6*23 octant degeneracy. Presence of the AP"°''™ ™^(z/^ —>■ u^) 
in first order in e indicates that the z/^ disappearance channel would play a role in lifting 
the sign-Amgj^ degeneracy if the measurement is done off the vacuum oscillation maximum. 

Now, we discuss the inverse problem, namely, whether the ^23 octant degeneracy can 
be resolved in the presence of the sign-Amg^ degeneracy. By noting that J^^*^* — J^^'^ = 
cos 26*2^1* Jr^*^* in leading order in cos 2^23, the difference between probabilities with the first 
and the second octant solutions can be given by 



AP 



1st 2nd/ 



= COS 2^23*^21 [sin^ 2^12 A21 + 4J^^^* (cos S sin 2A31 - 2 sin 5 sin^ A31)] 
+ API^* ^'^'^(z/^ ^ z/,; eer)Nsi + AP^^* '^''{u, ^ z/e; ee,)NSI, 



(74) 
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where 



AP 



1st 2nd , 



X 



l^e', ^erjNSI 



COS( 



-2v^Ci2Si2 COS 2^2'^ sin 26^^ 



f aL\ f Arni 



21 



A31 - -sin2A3i 



\eer\ sin0e^ siii^ A31 



(75) 



Apl«'2nd(^^ ^ i^e;^eM)iV5/ = 8 COS 2^23* 



X 



S23 Si3 



lee^l cos((5 + 0e^) ( — - — — + ^ siii 2A31 ) - l^e/xl sin((5 + 0e^) sin^ A31 



V A, 



31 



+ 



C12S12 / Am: 



21 



2v/2 V^"^^ 



2 J I |ee/x| cos0e/x 3A31 - - sin 26123 sin2A3i ) - l^e^l sin sin 26^21 sin A31 



(76) 



The first term of AP^*^* "^^^ in (!74|) . being composed only of the vacuum oscillation terms, 
is obviously invariant under the replacement normal inverted solutions. The remarkable 
feature of (!75|) and (!76ll is that they are also invariant under the replacement relation between 
different hierarchy solutions given in ( !66l) which is extended to include NSI phases. The 
disappearance channel does not play a role in the present discussion under the approximation 
taken in deriving (l65l) . because then AP^^* ^'^'^(z/^ — > v^) vanishes. Therefore, even in the 
presence of NSI, the resolution of the ^23 octant degeneracy can be carried out without 
worrying about the presence of the sign-Amg^^ degeneracy. The sign-Am^ and the 6'23 octant 
degeneracies decouple with each other even in the presence of NSI in matter perturbative 
regime. 



E. Non-Decoupling of Intrinsic degeneracy 



Now we discuss the intrinsic degeneracy for which the situation is somewhat different. 
First of all, this is the degeneracy which is somewhat different in nature. Unlike the case 
of the sign-Am|^ degeneracy, this degeneracy is known to be fragile to the spectrum anal- 
ysis; In many cases it can be resolved by including informations of energy dependence in 
the reconstructed events. An exam ple for this is the T2KK setting which receives an in- 
tense neutrino beam from J- PARC 0, S^l- It means that in this case there is no intrinsic 



degeneracy from the beginning. Nonetheless, anticipating possible circumstances in which 
spectrum informations are not available, and for completeness, we discuss below if resolving 
the intrinsic degeneracy decouple to lifting the other two degeneracies. We disregard the z/^ 
disappearance channel in this subsection because it does not appear to play a major role 
in resolving the intrinsic degeneracy. The discussions in this subsection are also meant to 
partly correct and append the ones given in Sec. Ill in 42 . 



1 . N on- Decoupling of Intrinsic degeneracy without NSI 



Let us first discuss the problem of decoupling with intrinsic degeneracy without NSI. 



In our perturbative approach AP 



12/ 



Uf,) arises only from the first order matter term 
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in (l62l) because the degenerate solutions in vacuum, by definition, gives the same vacuum 
oscillation probabilities. It reads 



where 



-4 



4e) 



1 2 A 1 

sm Aqi sm 2A-i-i 

A31 2 



(77) 



— sin 2^12 cot ^^23^21 cot A31 (^26^^r^ cos 5^^) + sin 2^12 cot 023 A21 cot A31) (78) 



Based on the result of AP^^ in (177|) we discuss possible decoupling of the sign-Am^ and the 
octant 623 degeneracies from the intrinsic one. 

We start from the sign-Am^ degeneracy. It can be readily seen that AP^^ is odd under 
interchange of the normal and the inverted hierarchy solutions as dictated in fl7ip . It means 
that AP^"^ (normal) — AP^"^ (inverted) = 2AP^^ . Clearly, the sign-Am^ degeneracy do not 
decouple from the intrinsic one. 

Now we turn to the octant ^23 degeneracy. From fl77|) . AP"^^(lst) — AP"'^^(2nd) reads 



AP^2(lst) - AP^2(2nd) = 8 

X cos 26^23 



ae) 



0«cos5«^, 



) sin 2^12 A21 cot A31 
1 + C23S23 



+ 



C23-S23(C23 + S23) siu 26*23 



1 2 A 1 

— sm A31 - -sm2A3i 

^31 ^ 



■ sin 2^12 A21 cot A31 



(79) 



where Oi^ and S23 etc. in ( !79l) are meant to be the ones in the first octant. It is small in the 
sense that it is proportional to cos 2^23 which vanishes in the limit of maximal ^23- But, this 
is the factor of kinematical origin which inevitably exists because the measure for breaking 
of the octant degeneracy has to vanish at 6*23 = 7r/4. Therefore, we conclude that there is 
no dynamical decoupling of the ^23 octant degeneracy from the intrinsic one. 

Now, we discuss the inverse problem, namely, whether the sign-Am^ and the ^23 octant 
degeneracies can be resolved independently of the intrinsic degeneracy. The measure for 
resolving the sign-Am^ degeneracy is given in (!72|) 



AP 



norm mv/ 



4A0^.^3 ( ^ 



norm mv 

1 



(2) 



A 



sin^ A 



31 



31 



- sin 2 A31 



(1) 



(80) 



where all the quantities in flHOl) is to be evaluated by using the normal hierarchy and intrinsic 
first solution. Clearly, the intrinsic degeneracy does not decouple from the sign-Am^ one. 
How about the 6*23 octant degeneracy? The appropriate measure for the question is given 

by 



AP 



1st 2nd/ 



1) - AP 



^^13 



1st 2nd ^2) =COS202\'*A21 



4J;" (^S + e^) cos5(^) + ^} sin2A3i - 2^^^Ae' sm' A,, 



^13 + ^13 



(81) 



where ^^3^ implies to insert the expression in (l69l) . Again there is no sign of the decoupling. 
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Nonetheless, there are some cases in which the decouphng with the intrinsic degeneracy 
still holds in a good approximation. For example, AP^^(lst) — AP^^(2nd) in (179!) and 
^pnorminv(^^) - AP^°™'°^ in §^ may be small numerically. It is the case at relatively 
short baseline L < 1000 km where it is further suppressed by The AP differences 
between the two 623 octant solutions are always suppressed by cos 26*23, and hence they may 
be small at 623 very close to the maximal. 

It is significant to observe that at the vacuum oscillation maxima, A31 = {2n + the 
decoupling is realized in all pairs of degeneracies. Therefore, if the experimental set up is 
near the vacuum oscillation maxima the decoupling with the intrinsic degeneracy perfectly 



holds. The identical two detector setting in T2KK [42, l52|, whose intermediate (far) detector 
is near the first (second) oscillation maximum provides a good example for such "accidental 
decoupling" . 



2. Decoupling and non- decoupling of Intrinsic degeneracy with NSI 



We concisely describe what happens in the decoupling between the intrinsic and the 
other two degeneracies when NSI is introduced. We explicitly discuss below the case with 
Est because the equations are slightly simpler, but we have verified that the same conclusion 
holds for the case with Eefj., and hence in the full system. 

^er type NSI gives rise to contribution to the difference of the probabilities with the first 
and the second solutions of intrinsic degeneracy of the following form 



X (^26^«cos(5« + ^^ 



C23S23 



/sin2A3i 

COS0er ' 



1 



V A. 



31 



sin 2 A31 I + sin (per sin^ A31 



(82) 



where use has been made of the relation ( ID 71) . Notice that the terms proportional to the 
solar Am2i do not contribute, and sin 5 terms cancel out owing to the relation (lD6p . 

We observe that AP^^(£:eT) are invariant under interchange between the normal and the 
inverted hierarchies, flTTl) . Therefore, NSI induced oscillation probability, by itself, fulfills 
the decoupling condition with the sign-Am^ degeneracy. 

The situation is different in relationship with the 6*23 octant degeneracy. With Eer one 
can derive the similar expression as fl79|) : 



(83) 



/\P^'^{Eer]l^t) - AP^2(^ee^. 2nd) = 4^2 COS 2^23 (^^^ l^erl sin2^i2A2iCot A, 

'sin2A3i 1 



COS( 



A 



31 



— - sin 2 A31 ) + sin sin^ A31 



Though the intrinsic degeneracy does not decouple with the 6*23 octant degeneracy, the 
suppression factor cos 26*23 (f§) kerl may be very small if baseline is relatively short and 6*23 
is near maximal, assuming the likely possibility that leerl is small. Again, the decoupling 
holds at the vacuum oscillation maxima. 

General conclusion in the last two subsections is that although the decoupling between 
the sign-Amg]^ and the 6*23 octant degeneracies holds, there is no decoupling between the 
intrinsic degeneracy and the other two types of degeneracies. The conclusion applies to the 
cases with and without NSI. 
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IX. CONCLUDING REMARKS 



In this paper, we have discussed various aspects of neutrino oscillation with NSI, the 
exactly hold properties as well as the properties best illuminated by a perturbative method. 
The former category includes the relation between the S matrix elements and the proba- 
bilities that arises due to an invariance of the Hamiltonian under the transformation flTB]) 
which involves 6^23 and the NSI elements Ea/s = e,fi,T). It allows us to connect the 

probabilities of various flavor conversion channels, which is powerful enough to strongly con- 
strain the way how various NSI elements Eajs enter into the oscillation probabilities. This 
category also includes the phase reduction theorem which guarantees reduction of number 
of CP violating phases when the solar Am^i is switched off. 

By taking the following three quantities, ■Sia, and the NSI elements Eoh, as small 

expansion parameters (which are collectively denoted as e) we have formulated a perturbative 
framework which we have dubbed as the "e perturbation theory" . Within this framework we 
have calculated the S matrix elements to order and derived the NSI second-order formula 
of the oscillation probability in all channels. It allows us to estimate size of the contribution 
of the particular NSI element Ea/s = e,fj.,T) to the particular oscillation probability 

P{^K ^Lj) {i^,uj = e,fi,T), as tabulated in Table [B To complete the table (and for other 
reasons) we have also calculated the oscillation probability in the z/g related channels to 
third order in e, which is given in Appendix O We have given a global overview of neutrino 
oscillation with NSI and hope that the table serves as a "handbook" for hunting NSI effects 
in neutrino propagation. 

Thanks to the NSI second-order formula we have discussed, for the first time, the way 
how the SI and the NSI parameters can be determined simultaneously. We found that 
measurement of all the relevant NSI and SI parameters is extremely demanding; While all 
the NSI elements in i/^ related sector can in principle be determined together with 613 and 
6, it requires I'e ^ ^fi, i^ti ^ ^e, ^ ^t, and their CP conjugate channels if we do it by the 
rate only measurement. We have also proven to the accuracy of that, if we restrict to the 
rate only analysis, all the NSI elements in z/^ — Ur sector cannot be determined even if we 
prepare Ur beam. 

Clearly, the right strategy is to pursue the appropriate experimental setup which enables 
us the spectrum analysis to determine several coefficients at the same time. The capability 
of spectrum analysis with good resolution would be a mandatory requirement for future 
facilities which aim at searching for effects of NSI at least as one of their objectives. To our 
knowledge, the leading candidate for such setup is the two-detector setup at L ~ 3000 km 
and L ~ 7000 km in neutrino factory with use of the golden channel [201] , Av hich are proven 
to be powerful in resolving the conventional parameter degeneracy [33|, l63j. In a previous 
paper, it was shown that the setting is also powerful in resolving the ^13-NSI (and probably 
the two-phase) confusion [31^. It must be stressed, however, that we still do not know if the 
setting is sufficiently powerful in determining all the SI and the NSI parameters. 

We have observed that the phenomenon of parameter degeneracy prevails in the system 
with NSI. Notably, it exists in an extended form of involving not only the SI but also the 
NSI parameters. In a concrete setting of six probabilities at monochromatic beam, we have 
uncovered a new type of degeneracy, the solar-atmospheric variable exchange degeneracy. 
To have a first grasp of the nature of the parameter degeneracy of more conventional type, 
we have discussed the matter perturbation theory of neutrino oscillation with NSI. We have 
found that the sign-Am|]^ and the 623 octant degeneracies are robust, and the analysis 
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indicates the way how the NSI parameters are involved into the new form of degeneracy. 
The decouphng between degeneracies, a sahent feature in the matter perturbative regime, 
is also revisited in an extended setting with NSI. 

In our investigation we have also noticed a new feature of neutrino oscillation in matter 
in the standard three-flavor oscillation without NSI, that is, the matter hesitation. It states 
that the matter effect comes in into the oscillation probability only at the second order in 
e. The property allows us to understand why it is so difficult to detect the matter effect in 
various long-baseline experiments, and explains why See is absent from the NSI second order 
formula. Notice that the property does not hold in the — Vr system with NSI. 

Of course, a number of cautions have to be made to correctly interpret our results; Many 
of our statements are based on the NSI second order formula which is reliable only if the 
assumptions made in formulating our perturbative treatment are correct. We do not deal 
with effects of NSI in production and detection of neutrinos. The program of complete 
determination of the NSI parameters mentioned above must be cooperated with search for 
NSI in production and detection processes. 

In this paper we confined the case of relatively small ^13 in accordance to our perturbative 
hypothesis in (fT9l) . What happens if ^13 is large enough so that not only 6*13 but also 5 are 
determined by the next generation reactor /accelerator [64, 65, 6^ and upgraded superbeam 
(67| experiments prior NSI search? Then, one might argue that the discussion of parameter 
determination would become much less complicated in this case. We argue that this is 
not quite correct. As we have seen in Sec. |Vl]the NSI and the SI parameters appear in the 
oscillation probability in a tightly coupled way. Hence, determination of the former with size 
of Ea^p ~ 10~^ requires simultaneous determination of the latter with accuracy of the similar 
order. Therefore, prior determination of ^13 and 5, unless extremely precise ones, would 
not alter the necessity of simultaneous determination of SI and NSI parameters. However, 
we note that knowing the neutrino mass hierarchy would greatly help by decreasing the 
ambiguities which arise from the degeneracy. 



APPENDIX A: S MATRIX ELEMENTS FOR NEUTRINO OSCILLATION 
WITH NSI 

Using the formalism described in Sec. |V] with the double-tilde basis fl29l) it is straight- 
forward to compute the S matrix elements for neutrino oscillations with NSI. Omitting 
calculations we just present the results of the 5* matrix elements: The notations used below 
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Am: 



'31 r , 



and the NSI elements are in the tilde-basis (fTSll . 



+ s^3(irAAL)e 



13 1 „ l*^ 



lAL 



1 - r^i 



) 



- 2si3Re(£ere''')rA 



1 - 



(^12 ^ee 



ci2^i2^ + eeMr{(^rAAL)e--^^ - (l - e"--^^) } 



l-TA 



l-TA 



(Al) 



-C23lCl2Si2^+£eM) (l - 6"^^^^^ 



) - S23 (Sise + rAEer) Z ( 



C23Si3£^.e - ^ (1 - e-*^^-^^) + (1 - e"*^^) } 



I-ta^ 



ir 



+ C23 C12S12 h ^eu 



rA 



2 ^ , ^ 



-ir^i AL 



+ .23 {sise-^' + e.r) [y^) 



'2 _ 2 n'^A 

rA 



2AL<! ( sl/^ + e 



I-ta 



'12 



rA 
''rA 



+ |c23£:*^ (Si3e + ^er) + S23^Mr (^Ci2Si2^ + ^e/.^ } 



X Ta 



:i-e 



-iAL\ 



1 - Ta 



(A2) 
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1 



Ser = S23 ( Cu-Su— + EeA {l - 6 '''^^^) - C23 (sigC + rASer) T-—- ( 
Ta J ^ 1 — Ta 



+ «23«i3£;.e-^{ (1 - e-^^-^^) + (1 - e"'^^) } 



+ C23Si3e '"^(e 



ee '~TTj 



l-TA 



+ S23 1 Cx2Sviy- + £e/. 



'12" 



rA 



-ir^ AL 



+ I (4 - ^L)^ - ^~ee + ) (1 - e-^^-^^) 



+ C23 (Sl3e + £er) ( Y 



ta 



-iAL 



l-TA 



2 rA 

rA 



^g-ir^AL _ ^-^AL^^ 



+ 
X Va 



:i-e 



-iAL\ 



l-TA 



l^g-ir^AL _ g-iAL^ 



(A3) 
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- 2C23S23 M^^.r)rA (1 - e-^^^) 



■iAL 



2 2 
^23*13 



1 - r A 



1 - r A ^ 



+ 2s23Sl3Re(£ere''^ 

+ 2C23S23S13 Re(£e,.e^')rA (l - e"^^^) 
+ 2c23S23Si3{ci2Si2COs5^ + Rc (eV^'') } (l - e"^^^^^) 

(rAAL)2 



rr>2 ^ , ~ N2 , 2 ~2 -iAL 
<-23V'"12„ ^^m) ^ °23^TT^ 



'23 



ci2«i2- + ee,?{{irA^L) - (1 - e--^^) | - (l - zAL - e^^^) 
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-iAL 
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ta 
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-itaAL g— iAL 



(A4) 
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Stt = S23{l - ^(c?2'^A + e^f,rA)AL^ + cl^ (1 - ierrVA^L) e 

+ 2C23S23 Rc(£V)^A (1 - e-^^^) 



2 2 
^23* 13 



1 - r A 



+ 2c23Sl3Re(£eTe 



i5N 



1 - r A ^ 



2C23S23S13 Re(£eMe'')^A (l - e"^^^) 

2c23S23Si3{ci2Si2 COS 5^ + Re {s.^e'^) } (1 - e-^'-^^^) 



+ *23 



-23 
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rA 



- 2c23S23Re|^Ci2Sl2^ +£e^^ (sisC + ^er) | 



X Ta 



1 - Ta 



(A5) 
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Sfj,T = -C23S23|l - iic^rA + £;u;,rA)AL| + C23S23 (1 - iSrrrA^L) e 
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l-VA 



— e 
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^—ivA^L —iAL 
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iAL 



Ta 



+ 
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ta 

1 _ e-^^^) - 
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TA 



TA 



C23S23 siae 



rA 

-i5 , 



l-TA 

2 



(< 



-iAL 



?2 -iAL 



iALe 



-iAL 
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^g-^r^AL _ g-iAL^ 



- C23S23|ci2Si2^ + ^^^^[irAAL - (l - e"^^^^^) } 



ta 

C23S23\e^.r?r\[iAL{l + e-*^^) - 2 (l - e" 



-irxAL 



)} 



The other 5" matrix elements are given by either the T-conjugate relations 

Sre{5,(f)ap) = SeriS, -(pa/s) , 

or by the CP-conjugate relations for antineutrino channels 



Seri-S, 



al3, —a), 
'af3, —a)- 



(A6) 



(AT) 



(A8) 
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APPENDIX B: NSI SECOND-ORDER PROBABILITY FORMULAS 



In this Appendix we give the exphcit expressions of the oscillation probabilities to second 
order in e in all channels, except for those which can be readily obtained by the extended 
transformation flTBl) . 

1. Oscillation probability in the i/g-related sector 

We present here the explicit forms of ^(t'e ^ ^'e) and P{i'e ^fi) for completeness and 
possible convenience of the readers considering importance of the appearance channels. 
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^^21 ^/ i5^"^31 * 
IC12S12 h C23£e,i - S23£^eT)(Sl3e ^ V 8236^^ + C23£e 



a 



a 



aL 



X 



Aml^ 



AmoiL Ami, - 
sm — - sm — -— — sin -— 

AE AE AE 



-L.(B2) 



PiVf, Vt) can be obtained from Piv^^ u^) by the transformation f|T6|) . Or, the simpler 
way of remembering the operation is to do transformation C23 —S23 and S23 — >■ C23 in 
P(z/e — * z/^), but undoing any transformation in the generalized atmospheric and the solar 
variables defined in fl55]) . the pieces bracketed in the real and imaginary parts in flB2p . (See 
also (!B3|1 . 



2. Oscillation probability in the i/^ — Ur sector 

For compact expressions of the oscillation probabilities in the — u^- sector, we define 
the simplified notations which involve e's in the z/^ — Ur sector as well as ^ee-^^ Together 



19 readers who want to see the fully explicit expressions of all the oscillation probabilities, we refer the 
first arXiv version of this paper 68 1. 
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with the ones already defined in Sec. IVI Al they follows: 

C12S12 h C23€ef, - S23£eT 6 = ^ 6 S 

a J 

= C23S23(£^^ — Err) + C^^^^r — sls^^tr = \^\^^'^ i 

Si = (C23 - S23)(£rT - ^^mm) + 2C23S23(£^r - C?2— 

'^2 = {S^i^ — See) + C^^i^rr — ^fifj.) + C23'523(^/xt + ^lir) ~^ ^12 

'^3 = {Sfifj. — See) + sl^i^TT — Sfifj.) — C23S23{ef^r + sI,t) + {^12 ^ ^12) (^3) 

Notice that Si, S2, and ^3 are not independent, Si = S2 — S3. We also note that 6±, S, 
and £^ are complex numbers while the others are real. 

To present the oscillation probabilities in the — v^- sector, we start by recapitulating 
the decomposition formula (146!) in Sec. IVIBt 

P(z/a vp] Ee^, Ser, £ £fiT, £tt) = P{i^a ^^p] 2 fiavor in vacuum) 

+ P{^a tyi3]£efj.,£eT) 

+ P{i^a ^ i^l3]£f,f„£f^T,£TT) (B4) 

where a and jS denote one of /i and r. The first term in (IB4p has a form that it appears in 
the two fiavor oscillation in vacuum: 

P(z/^ v^j_\2 fiavor in vacuum) = P(z/^ ^ z/,-; 2 fiavor in vacuum) 

,2 „2 2 Am|^L 



1 — 4cnoSno sin 



23023 0111 , 

P{^fi ~^ ^t] 2 fiavor in vacuum) = 4C23S23 sin^ — — . (B5) 
We have shown in Sec. IVI B I that the third term in (lB4p in the ^ u^, v^- ^ u^-, and 
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z/t- channels are given by the single equation' 
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The second term in (lB4p is given in z/^ — z/^ channel as 
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The subscript ± in this and the following equations denotes the normal and the inverted 
mass hierarchies, which corresponds to the positive and negative values of Am|]^. Notice 
again that P(z/t- — > z/,-; ^e/^, ^er) can be obtained from P(z/^ —>■ i'fj,;eefj,,eeT) by the extended 
transformation (|T6i) . or by the operation described at the end of the previous subsection. 



To second order in e the sensitivity to e^^ and Err is through the form e^^ — err; and hence no sensitivity to 
the individual e's. GeneraUy, the diagonal e's appear in a form of difference in the oscillation probabilities 
as one can observe in the third-order formula given in Appendix [C] It must be the case because the 
over-all phase is an unobservable. 
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Finally, the second term in the oscillation probability in the Vr channel is given by 
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APPENDIX C: NSI THIRD-ORDER FORMULA 

We present here the third-order formula for the oscillation probability with NSI. Though 
utility of such lengthy formula may be subject to doubt we can offer at least three arguments 
to justify the presentation of the formula in this Appendix. Firstly, the third-order formula 
is needed to complete Table [B Secondly, by turning off all the NSI elements one can obtain 
the SI third-order formula for the oscillation probabilities with SI only, which is valid to 
order e^. To our knowledge, such formula has never been derived in the literature. Utility 
of the SI third-order formula for theoretical analysis may be obvious if the sensitivity to 
the oscillation probability reaches to the level of ~ 10~^, which is smaller than terms of 
order e^. In fact, it appears to be the case in some of the future facilities according to the 
analysis in j69|. Thirdly, once the sensitivity to the oscillation probability comes down to 
~ 10~^, a complete treatment of neutrino oscillation probability must include NSI elements 
up to the same order as ^13 and Am^i/Amlj^, as far as our ansatz (fT9l) in formulating the 
e perturbation theory is correct. Thus, we believe that the NSI third-order formula has a 
good chance to be useful. 

In presenting the third-order probability formula we restrict ourselves to the z/g related 
channel, and only present P(t'e ^ii) here because from which P{Vf, Ut) can be obtained 
by the extended transformation (ITB]) . Then, P(z/e — u^) can be readily calculated by using 
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the unitarity relation. The NSI third-order formula for P(i/e — ^ J^ij) reads 
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APPENDIX D: INTRINSIC DEGENERACY IN VACUUM 



We re-examine the problem of intrinsic degeneracy in vacuum. For simplicity, we focus 
on the channel z/^ —>■ z/g. We use a simplified notation S13 = s below. The neutrino and 
anti-neutrino oscillation probabilities in vacuum are given by 

P{j^ii J^c) = Xs^ + (Fc cos 5-Ys sin 5) s + Pq 

^ yc) = Xs^ + (Fc cos 5 + Ys sin 5) s + Pq (D1) 

Am-? L 

where X, F's, etc. are defined with simplified symbol Ajj = as 

X = 4s23sin^A3i, 

Yc = sin 26'i2 sin 2^23^21 sin 2 A31, 

Ys = 2sin2^i2sin2^23A2isin2 A31, 

Pq = sin^ 20124^21- (D2) 
Let us denote two set of intrinsic degenerate solutions as (si, Si) and (s2, 62). They satisfy 

P - Pq = Xsl + (Fc cos 61 - Ys sin 61) Si 

P-Pq = Xsl + (Fc cos §2 - Ys sin ^2) S2 (D3) 

and 

p - = Xsl + (Yc cos 61 + Ys sin 61) Si 

P-Pq = Xsl + (F, cos 62 + Ys sin 62) S2 (D4) 
By subtracting two equations in (IDSp and (]D4p respectively, we obtain 

X(s^ — S2) + Fc(si cos 5i — S2 cos 82) — Fs(si sin (5i — S2 sin ^2) = 0, 

— S2) + Fc(si cos (5i — S2 cos(52) + Fs(si sin 5i — S2 sin ^2) = 0. (D5) 
They further simplifies to 

si sin 5i — S2 sin (^2 = 0, (D6) 
X{s\ - si) + Fc(si cos 5i - S2 cos 82) = 0. (D7) 

Equation (]D6p can be solved as 



S2 cos ^2 = ± Y ■§2 — sin 61 (D8) 

which can be inserted to (]D7|) to yield the (formally quartic but actually) quadratic equation 
for S2. Now, the issue here is to choose the correct sign in (IDSp . One can show that by using 
( ]D7p if Fc > (Fc < 0), minus (plus) sign has to be chosen. 

These equations can be easily solved for (s2, ^2) for given values of (si, 61) as inputs: 



F\ . /F^ ' 



■S2 = \lsl + 2 Si COS^i + y ^ 

sin (^2 = — sin 61 

cos52 = =F— ( si cos5i + ^ ) (D9) 
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where the sign =p for cos(52 are for Yc = ±|^c|, and S2 in the solution of S is meant to be the 
S2 solution given in the first line in (lD9p . By using 



Y 

= sin 2^12 cot ^23A2i cot A31 (DIO) 



S2 can be written as 



■S2 = Y + 2 sin2^i2 cot ^^23A2i cot A31S1 cos5i + (sin26'i2 cot 6'23A2i cot A31) (Dll) 
Similarly, cos 6 is given as 

cos(52 = =F — (si cos (5i + sin 2^12 cot ^23A2i cot A31) (D12) 
By further expanding (1D9P by ^, assuming it small, the Burguet-Castell et al. solution 



33l . l57l | is reproduced; 



S2 ^ si + cos5i (D13) 
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